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D , Abstract. It is well known that systems with long-range interactions may exhibit 

^ ■ different phase diagrams when studied within two different ensembles. In many of the 

previously studied examples of ensemble inequivalence, the phase diagrams differ only 
^ I when the transition in one of the ensembles is first order. By contrast, in a recent study 

■ of a generalized ABC model, the canonical and grand-canonical ensembles of the model 

were shown to differ even when they both exhibit a continuous transition. Here we 
show that the order of the transition where ensemble inequivalence may occur is related 
to the symmetry properties of the order parameter associated with the transition. 
' ^ I This is done by analyzing the Landau expansion of a generic model with long-range 

Q ■ interactions. The conclusions drawn from the generic analysis are demonstrated for 

. the ABC model by explicit calculation of its Landau expansion. 
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diffusive systems (Theory) 
1. Introduction 



CN| ! Much work has been devoted in recent years to the understanding of the properties of 

systems with long-range interactions. These are systems where the two-body potential 
decays at large distance, R, as 1/R'^~", with < a < ci in ci dimensions. Interactions of 
^ ! this type are found in various systems such as self-gravitating models [1, 2], non- neutral 

5^1 plasmas [3], dipolar ferroelectrics [4] and others. 

The algebraic decay of the two-body potential results in an ill-defined 
thermodynamic limit since the energy is super-extensive {E ~ y^+f^/d^ where V is the 
volume of the system). A thermodynamic limit may be restored by applying the Kac 
prescription [5], by which the potential strength is rescaled by an appropriate volume- 
dependent factor. However, even within this scheme the energy remains nonadditive, in 
a sense that the energy of a system composed of two subsystems is not necessarily equal 
to the sum of the energies of the two. Nonadditivity results in some distinct properties 
that are not found in the more commonly studied short-range interacting systems. These 
properties include, for example, negative specific heat in micro-canonical ensembles (and 
similarly negative compressibility in canonical ensembles) and inequivalence of statistical 
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ensembles, which is the focus of this paper [6-19]. Long-range interactions may also 
result in dynamical effects such as the existence of quasi-stationary states, which are 
long-lived states different from the usual Gibbs-Boltzmann one [20, 21], and anomalous 
relaxation time to equilibrium. Recent reviews of this topic are found in [22-26]. 

In many studies of ensemble inequivalence in models with long-range interactions, 
it was found that whenever the two ensembles (either micro-canonical/canonical or 
canonical/grand-canonical) yield a continuous transition, the transition occurs at the 
same critical point. On the other hand, the phase diagrams of the two ensembles 
differ from each other when the phase transition becomes first order in one of the 
ensembles. Such behaviour has been observed in studies of numerous systems such as 
the spin-1 Blume-Emery-Griffiths (BEG) model [13, 14, 27, 28], the ABC model [17, 18], 
the zero range process [16], two-dimensional vortices [29], random field and spin glass 
systems [30-32] and others. However, general considerations have shown that different 
ensembles may, in principle, exhibit different critical points [33]. In fact, a recent study 
of a generalized ABC model revealed a case where the canonical and grand-canonical 
ensembles exhibit a second order phase transition at different critical points [19]. 

In this paper we consider a general framework within which one can relate the 
order of the phase transition where ensemble inequivalence may be observed and the 
symmetry properties of the underlying order parameter of the transition. To this end, 
we consider a generic model which undergoes a phase transition governed by some order 
parameter, m. This parameter could, for instance, be the average magnetization in the 
case of a magnetic transition such as in the Ising model. We consider the model within a 
'higher' ensemble, where a certain thermodynamic variable, denoted by q, is allowed to 
fluctuate, and within a 'lower' ensemble, where q is kept at a flxed value. For instance, 
in the case where q is the energy the two ensembles would correspond to the canonical 
and micro-canonical ensembles, whereas in the case where q is the particle density they 
would correspond to the grand-canonical and canonical ensembles. 

A natural framework within which phase diagrams of systems with long-range 
interactions can be analyzed is the Landau theory. In this approach one considers 
the partition sum over the microscopic configurations of a system, C, that correspond 
to some value of the thermodynamic variables, taken here to be m and q, given by 

eM-PVf{m,q))= J] exp(-/3^(C)). (1) 

C : m{C)=m,q{C)=q 

Here E is the energy of a configuration, V is the volume of the system and /(m, q) is 
the large deviation function. Landau theory is a phenomenological approach whereby 
f{m, q) is assumed to be an analytic function even at criticality. This allows one to study 
the phase diagram of the model by expanding /(m, q) in powers of m and q around 
its minima. In systems with long-range interactions, the Landau theory is expected 
to be exact to leading order in 1/V. This conjecture, whereby f{m,q) is analytic, 
is of course true for models with infinite-range mean-field type interactions, and has 
been proven rigourously also for spin models with l/R'^^'^ interactions for < a < o? 
under certain conditions [34, 35]. We can therefore study the behaviour of long-range 



3 



interacting models in a rather general context by considering a generic Landau expansion 
of f{m, q). Since the symmetry properties of m and q determine which terms may appear 
in the Landau expansion, this approach allows us to study the connection between the 
symmetry of the model and the nature of its phase transitions. 

In the present paper we show how the symmetry properties of m and q determine 
whether or not the two ensembles exhibit the same critical points. For simplicity we 
study here the case where m and q are scalars and where f{m,q) is invariant under a 
simple transformation of m, of the form m — > — m. Our results can be straightforwardly 
generalized to more elaborate cases, where for instance m is a vector or a complex 
number, as demonstrated below in our study of the ABC model. In the 'lower' ensemble, 
where q is fixed, the Landau expansion of /(m, q) around m = is given by 

/(m, q) = a(q) + c(q)rn' + e(q)m^ + O(m^). (2) 

The critical point of the 'lower' ensemble is found where the coefficient of changes 
sign, and it is thus given by the equation c{q) — 0. In the 'higher' ensemble, fiuctuations 
in q around its average value g*, denoted by 5q — q — q*, may or may not affect the 
coefficient of depending on the symmetry of the model. In order to study this effect 
we expand the coefficients in (2) in powers of 5q. Two types of symmetries are considered 
in this expansion: 

(i) f{fn,q) = f{—ni,q) : In this case the dominant coupling term between m and Sq 
in the Landau expansion is of the form m^5q. As a result, near the critical point of 
the 'higher' ensemble 5q scales as 5q ~ m^, implying that 5q affects the coefficient 
of the rn^ term in the Landau expansion and higher order terms. Since the critical 
line is determined by the coefficient of m?, the 'lower' and the 'higher' ensembles 
yield the same critical line. 

(n) /(m, q) = f{—m, —q) : In this case we study the stability of phase where q* — for 
which this symmetry is relevant. Here f{m, q) may involve a bilinear term of the 
form mSq, which was not considered in (2). As a result, near the critical point of 
the 'higher' ensemble Sq scales as Sq m. This scahng implies that Sq affects the 
coefficient of the term in the Landau expansion, causing the 'higher' ensemble 
to exhibit a critical line which is different than that of the 'lower' ensemble. 

The first symmetry, f{m,q) — f{—'m,q), corresponds to the more commonly studied 
case, where for instance m is the average magnetization and q is either the energy of 
the overall particle density. This is because both the energy and the particle density are 
invariant under any symmetry transformation of the model. 

The conclusions drawn from the above generic analysis are demonstrated explicitly 
for the ABC model [36-38]. This is a three-species exclusion model defined on a discrete 
one-dimensional lattice. The model exhibits a phase transition between a homogenous 
state and a phase-separated state, which is composed of three coexisting regions each 
predominantly occupied by one of the species. Here we study the ABC model on an 
interval, with closed boundary conditions. In this case, although the model evolves 



4 



by local dynamics, its steady state is given by the Gibbs-Boltzmann distribution of a 
system with long-range interactions [38]. 

The original ABC model evolves under particle-conserving dynamics where the 
overall density of each species, denoted by r^i, and tq, is fixed. Previous studies have 
introduced two kinds of particle-nonconserving generalizations of the ABC dynamics, 
corresponding to two distinct grand-canonical ensembles. In the first generalization, 
referred to as Model 1, the overall density of the three species, r = r^i -|- -|- r^, is 
allowed to fiuctuate while the differences between their densities, ta — tb and tb — rc, 
are fixed. In the second generalization, referred to as Model 2, the differences, — tb 
and rs — rc are allowed to fluctuate while maintaining the overall particle density, 
r, constant. Here we show that the two grand-canonical ABC models obey different 
symmetries, corresponding to the two different generic models mentioned above. In 
the ABC model the symmetry transformation corresponding to m —m is a cyclic 
exchange of the labels of the three species, A ^ B, B ^ C and C ^ A. Since the 
fluctuating fleld of Model 1, r — ta + rs + rc, is invariant under this transformation. 
Model 1 and the canonical ABC model exhibit the same critical fine, as shown previously 
in [17, 18]. On the other hand, since the fluctuating flelds of Model 2, rA—rB and rB—rc, 
are not invariant under the label-exchange. Model 2 and the canonical ABC model 
exhibit different critical lines, as shown previously in [19]. By carrying out the Landau 
expansion of the two grand-canonical ABC models explicitly, we derive their critical 
points and demonstrate the two models are related to the generic models mentioned 
above. In addition, the Landau expansion of Model 2 provides an alternative derivation 
of its critical point, computed previously in [19] by deriving the exact expression of the 
steady-state density proflle of Model 2. 

The paper is organized as follows: In section 2 we study ensemble inequivalence 
in the context of two generic models with long-range interactions, each exhibiting a 
different symmetry, as discussed above. The results of section 2 are demonstrated on 
the ABC model, which is presented in section 3.1. We study the two grand- canonical 
generalizations of the ABC model, Model 1 and Model 2, by explicitly computing 
their Landau expansion in sections 3.2 and 3.3, respectively. Concluding remarks are 
presented in section 4 

2. Ensemble inequivalence in a generic Landau model 

2.1. Symmetry of the type f{m,q) — f{—m, q) 

In this section we consider a generic model whose Landau free energy is invariant under 
the transformation {m,q) — > (— m, q), such as in the case where m and q correspond to 
the overall magnetization and the energy of the system, respectively. We study the phase 
diagram of the generic model within a lower ensemble, where q is fixed, and within a 
higher ensemble where q is allowed to fiuctuate, and show that the two ensembles exhibit 
the same critical line. 
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Lower ensemble: In the lower ensemble, where q is fixed, and close to the transition 
temperature where m is small, the Landau free energy density can be expanded in 
powers of m as 

/(m, q) = ao{q) + Co(g)m^ + eo(g)m^ + 0{mP). (3) 

The above coefficients are labeled with subscript in order to be consistent with their 
expansion in small fiuctuations in q, carried out below in the study of the higher 
ensemble. The free energy of the lower ensemble is obtained by minimizing the Landau 
free energy with respect to m, 

F(r, ?,...)= min/(m, q) = /(m*, g), (4) 

m 

where m* is the average magnetization for a given value of q. Here T is the temperature 
and the dots denote other possible parameters of the system. In order simphfy the 
notation we omit here and below the dependence of the coefficients of the Landau 
expansion on T and those additional parameters. 

The phase diagram of the lower ensemble can be derived from (3) as follows: At 
Co = and for cq > 0, the model undergoes a second order transition between a 
disordered phase, for which Cq > and m* = 0, and an ordered phase, where Cq < and 
m* 7^ 0. This transition fine is denoted by the thin sohd line in figure la. The figure 
shows a sketch of the phase diagram of the model in the (co. A) plane, where A is the 
conjugate field of q, given by 

X{q) = df{m*,q)/dq. (5) 

On the second order transition line X{q) changes smoothly. This transition line 
terminates at tricritical point (★), where cq = eo = 0. For cq < 0, the critical line 
is preempted by a first order phase transition, taking place at Cq > 0, whose value is 
determined by the higher order terms in the Landau expansion. The first order phase 
transition is denoted by the vertical lines in figure la, which mark the discontinuity in 
A(g), observed as the system switches between a phase where m* = 0, and a phase with 
a nonvanishing value of m*. 

Higher ensemble: In the higher ensemble, where q is allowed to fluctuate, the Landau 
free energy density is given by the Legendre transform of f{m, q), defined as 

g{m, A) = min [/(m, q) — Xq] . (6) 

The free energy of the higher ensemble is obtained by minimizing /(m, q) — Xq with 
respect to q and m while keeping A fixed, 

G(T, A, . . .) - min (/(m, q) - Xq) = /(m*, g*) - Xq\ (7) 

m,q 

Close to critical fine we consider small fiuctuations around q*, denoted hy q — q* + 5q, 
which yields the following expansion, 

f{m,q)-Xq= [a^iq*) - Xq*] + a^iq^i^qf (8) 
+ [coiq') + ci{q*')dq + C2iq')i6qf] + eo{q*)m^ + O(m^), 
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Figure 1. Sketch of the phase diagram of the lower (a) and higher (b) ensemble of the 
generic model where /(m, q) = f{—m, q). Both diagrams show the same second order 
transition line (thin solid line) which terminates at a different tricritical point (★) in 
(a) and in (b), where the transition becomes first order. The first order transition is 
denoted by a thick solid line in (b) and by a discontinuity in A in (a) (vertical lines). 



where an{q) = for a = a,c, e. Here, we omitted the hnear term in 5q since (8) is 

minimal at q = q* and m = 0, which imphes that 

ai(q*)^Xq\ (9) 

In order to obtain the Landau expansion of the higher ensemble we minimize (8) with 
respect to Sq, yielding 

Ci{q* 



Sq^ 



-m' + 0(m^). 



(10) 



2a2(g*) 

The scaling of the fluctuations in q, Sq ~ m^, is a direct consequence of the symmetry 
of f{m, q), by which the lowest order coupling term between m and Sq in /(m, q) is of 
the form irP'Sq. This scaling implies that the fluctuations in q affect only the coefficient 
of rn^ and higher order terms in g{m^ A). The Landau expansion of the higher ensemble 
is obtained by inserting (10) into (8) and the result into (6), yielding 



g{m, A) = Co(g*)m + 



eo(?*) - 



4a2(g* 



m^ + 0(m*'). 



(11) 



Here the dependence on A is expressed through ^^(A), given by (9). The higher ensemble 
therefore exhibits the same critical line as the lower ensemble, at Cq = 0. Here, however, 
the second order transition line terminates at a different tricritical point, given by 
eo = Ci/4a2. A sketch of the resulting phase diagram is shown in figure lb. The 
plot shows a second order transition line (thin solid line) which turns into a first order 
(thick solid line) at a tricritical point (ik) , which is different than the tricritical point of 
the lower ensemble. 

Figure 1 corresponds to the case where the coefficient of the term in /(m, q) 
and g{'m, A) is positive, and thus the second and the first order transition lines join at 
a tricritical point in both ensembles. When the coefficient of is negative the two 
lines join at a critical-end-point, yielding a different phase diagram which will not be 
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discussed here. An example of a critical-end-point was found in a variant of Model 1 in 
[18]. 

2.2. Symmetry of the type f{m, q) = f{—m, —q) 

In this section we consider a case where the generic Landau free energy has a symmetry 
of the form f{m,q) = f{—m, —q). A simple example of this symmetry is provided by 
the anisotropic XY model with infinite-range interactions, as discussed at the end of this 
section. This type of symmetry is also relevant for Model 2 studied in section 3.3. As in 
the previous section we study the generic model within a lower and a higher ensemble, 
which in this case yield different critical points. 

Lower ensemble: In the lower ensemble, where q is fixed, and close to the transition 
temperature, we consider a general Landau expansion of the form 

/(m, q) = ao{q) + h{q)m + Co(g)m^ + do{q)rin^ + eo(g)m^ + O(m^). (12) 

In order to satisfy the symmetry, f{m,q) = f{—m,—q), the coefficients of the even 
powers of m must obey 

a-oiq) = a,o{-q) Co{q) = Co(-g) eo{q) = eo(-g), (13) 
while coefficients of the odd powers of m obey 

bo{q) = -boi-q), o?o(g) = -do{-q). (14) 

We study the model close to the q — fine, where (13) and (14) play a role in determining 
which terms appear in the Landau expansion of the higher ensemble, as discussed below. 
For g = all the coefficients of the odd powers of m in (12) vanish, and the model 
exhibit a critical point at Co(0) = (note that Cq depends on T and possibly other 
parameters that are not mentioned here explicitly). The order of the transition at this 
point depends on the sign of eo(0). The phase diagram of the lower ensemble for the 
case of eo(0) > 0, where the transition is of second order, is sketched in figure 2a. 
It is easy to see from (12)-(14) that above the transition, for cq > 0, the q — fine 
corresponds to X{q) — df{m*, q)/dq — 0. Below the critical point (•), namely for cq < 0, 
the lower ensemble undergoes a first order transition at g = 0, in which the conjugate 
field of q exhibits a discontinuity as the system switches between a phase with a positive 
nonvanishing m* and a phase with a negative nonvanishing m* (vertical lines) . 

Higher ensemble: In the higher ensemble, where q is allowed to ffuctuate, we follow 
the same lines of derivation presented in section 2.1. We consider the expansion of 
/(m, q) — Xq in powers of m and Sq around q* = 0, which is given by 

/(m, 5q) = ao(0) + a2(0)(5g)2 + bi(0)m5q + Co(0)m2 + O(m^). (15) 

By minimizing / with respect to Sq we obtain to lowest order 
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(a) Lower ensemble 



(b) Higher ensemble 



6 ^c^ 

Figure 2. Sketch of the phase diagram of the lower (a) and higher (b) ensemble of the 
generic model where f{m,q) = f{—m,—q). Here the two ensembles display different 
critical points (•) at cq = in (a) and at cq = h\/Aa2 in (b). For lower values of cq we 
find a first order transition, denoted in (b) by a solid line and in (a) by a discontinuity 
in A (vertical lines). 



This hnear scahng, 5q ~ m, is a direct consequence of the symmetry of f{m,q), which 
aUows for a bihnear coupling term of the form m6q. As a resuh of this scahng, the 
fluctuations in q affect the coefficient of in g{m,X). Inserting (16) into (15) and the 
resuh into (6) yields 

m 



g{m,\) = ao(0) + 



co(0) - 



m' + 0{rrf). (17) 



4a2(0). 

We therefore conclude that the higher ensemble has a critical point at cq = 6^/4a2, 
which is different from the one found in the lower ensemble, given by cq = 0. The order 
of the transition at the critical point depends on the higher order terms in (17), which 
are not considered here explicitly. The resulting phase diagram of the higher ensemble 
is sketched in figure 2b for the case where the coefficient of in (17) is positive at the 
critical point. The figure shows the critical point (•), below which the higher ensemble 
undergoes a first order phase transition at A = 0. Here of course, since A is the control 
parameter, it does not exhibit a discontinuity. 

A simple example for a model that obeys the symmetry discussed above is the 
anisotropic XY model with infinite-range interactions. The model is defined on a one- 
dimensional lattice of N sites. The equihbrium state of the model is a Gibbs-Boltzmann 
distribution corresponding to the following Hamiltonian, 

\i=\ J \i=l / \i=l j=l / 

Here = ±1 are the x and y components of the spin at site i, respectively, 

and Jxx- ^ Jyy > 0- By defining the average magnetization in each direction as 
ma = ^ Sa,i for a = x,y we obtain a Landau free energy density of the form 

/(m^, my) = - ^ {Jxxml + Jxym^^my + Jyyml) (19) 
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l + nia-, . 1 + m„ , , 1 - m« , ,1 - rUc 

n n ) n "^i n 



a=x,y 

where we took ks — 1 in order to simplify the notations. 

We study the model within a lower ensemble, where rUy is fixed, and within a 
higher ensemble where rriy is allowed to fiuctuate. From the symmetry of (19), whereby 
f{mx,my) = /(— m^, —rriy), we expect to find a phase diagram similar to figure 2, with 
rrix and rriy playing the role of m and q, respectively. In the lower ensemble we follow 
the lines of derivation presented for the generic model and obtain a critical point at 
rUy — and 

T/-^ = Jxx. (20) 
In the higher ensemble, the Landau free energy is given by 

g{m^, Hy) = min [/(m^,, rUy) - HyVfiy] , (21) 

my 

where Hy is the magnetic field in the y direction, corresponding to field A of the generic 
model. By analyzing the Landau expansion of the model around rUx — rUy — and 
Hy — 0, we find that a specific linear combination of rrix and rUy becomes unstable at 



rpHE 

2 



2 i^Jxx + Jyy + {Jxx Jyy)^ + Jxy^ ■ (22) 

The resulting phase diagram is similar to the one sketched in figure 2. As expected, due 
to the symmetry of the model, whereby f{mrc,my) = /{—rrix, —my), the lower and the 
higher ensembles exhibit different second order transition points given by (20) and (22), 
respectively. 

Below the critical point, the XY model undergoes a first order phase transition. 
The nature of this transition can be understood by plotting the magnetic field in 
the y direction in the lower ensemble, given by Hy = df{m*,my)/dmy, where m* 
the solution of df{mx,my)/dmx — 0. In figure 3a, where T > T^^ , we find a one- 
to-one correspondence between Hy and rUy which implies that the two ensembles are 
equivalent. For lower temperatures, in figure 3b where T/^^ < T < T^^, since Hy 
becomes nonmonotonous the higher ensemble undergoes a first order phase transition 
at Hy = 0, signified by a discontinuity in m*y{Hy). On the other hand, the lower ensemble 
does not undergo any phase transition in figure 3b. For T < T/"^, shown in figure 3c, 
the lower ensemble also exhibits a first order phase transition as we vary rriy through 
rUy = 0, yielding a discontinuity in Hy. 



3. ABC model on an interval 



The general discussion presented in the previous section may be applied to the ABC 
model. This is a three-species exclusion model studied here on an interval, where 
its dynamics obeys detailed balance with respect to a Hamiltonian with long-range 
interactions. In section 3.1 we review the phase diagram of the canonical ABC model, 
where the overall densities of each of the three species are conserved. We then study two 
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Figure 3. The magnetic field in the y direction in the XY model, Hy = 

df{m*,my)/dmy, as a function rriy for J^^ = ^,Jyv = 0-5 and Jxy = 0.1. For this 
choice of parameters T^^^ = 1 and T^^ « 1.06. In (a), where T > T^^, the one- 
to-one correspondence between Hy and my implies that the ensembles are equivalent. 
In (b), where T^^ < T < T^^, since Hy becomes nonmonotonous only the higher 
ensemble undergoes a first order phase transition at Hy = 0, where it exhibits a 
discontinuity in niy. For T < T^^ , shown in (c), the lower ensemble also undergoes a 
first order phase transition at rUy = 0, where it exhibits a discontinuity in Hy. 

distinct grand-canonical generalizations of the model in sections 3.2 and 3.3, referred to 
as Model 1 and Model 2 , respectively, where different conservation rules are relaxed. 
By explicitly calculating their Landau free energy. Model 1 and Model 2 are shown to 
correspond to the generic models presented in sections 2.1 and 2.2, respectively. 

3.1. ABC model on an interval: canonical ensemble 

The ABC model is a one-dimensional driven exclusion model which has been studied 
extensively in recent years, mainly as an example of a model which, when studied on 
a ring geometry, exhibits long-range order and a phase-separation transition out of 
equilibrium [36, 37, 39]. Here we consider the model on an interval, where its dynamics 
satisfies detailed balance [38]. Each of the L lattice sites can be occupied by one particle 
(either A,B or C) or remain empty. The microstate is thus denoted by ^ = {Cij^li; 
where i is the index of the site and Q = A,B,C or 0. The numbers of A, B and C 
particles are denoted by iV^, Nb and Nc, respectively, with Na+Nb+Nc — N < L. We 
also denote the average density of each species hy ra — N^/ L ior a — A, B,C and the 
overall particle density by r = rA+TB+rc < 1- The model evolves by random sequential 
updates whereby particles on neighbouring sites are exchanged with the following rates, 

ABi^BA, BC^CB, CA^AC. (23) 
111 

The vacancies are inert and exchange symmetrically with the other particles with the 

following rates, 

AoioA, SoioS, Coioa (24) 
111 
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For q = 1 the model relaxes to a homogeneous steady state where the particles are 
randomly distributed on the lattice, whereas for q ^ 1 the particles phase separate into 
three domains in the limit L — )■ oo. The vacancies remain homogenously distributed 
in the system for any value of q. Without loss of generality we consider here g < 1 
where the domains are arranged clockwise as AAA . . . BBB . . . CCC, with vacancies 
distributed homogenously on the lattice. A unique feature of the ABC model on an 
interval is that although its dynamics is local, the model obeys detailed balance with 
respect to an effective Hamiltonian with long-range interactions, given by 

L-l L~i ^ 

^ (^) = E E (^^^^+^- + + ^^^^+^) - 6^'- ^^^^ 

i=i j=i 

The projection operators in the Hamiltonian are defined as 

/>.N J 1 for = a 
"'(^) = \0 else ' 

for a = A,B and C. The fact that a particle on site i interacts equivalently with all 
the particles to its right (on sites i + j with j = 1, . . . ,L — i) implies that "H (C) is a 
long-range Hamiltonian. For convenience, we add a constant term to the Hamiltonian, 
— |A^^, which sets the average energy of the disordered phase to zero. The fact that 
this terms scales super-linearly with the size of the system is another evidence for the 
long-range character of the ABC model. On a ring geometry, not studied here, the 
model obeys detailed balance with respect to (25) only for Na — Nb — Nq, whereas on 
an interval this is true for arbitrary Na, Nb and Nc- 

The ABC model is often studied in the limit of weak asymmetry, where q approaches 
1 in the thermodynamic limit as q = e~^^^. Here the parameter /5 serves as the inverse 
temperature of the model. For equal densities, A'^ = Nb = Nc, the model undergoes a 
second order transition in the L — )■ oo limit between a disordered phase and an ordered 
phase, with three macroscopic domains each predominantly occupied by one of the 
species. The transition point, given by 

(3, = 27iV3/r, (27) 

has been derived for r = N/L = 1 in [38, 39] and for r < 1 in [17]. For nonequal densities, 
the model relaxes to an ordered steady-state profile for all values of /3, and hence no 
phase transition takes place [38]. The properties of the weakly- asymmetric ABC model 
close to its transition point has been the focus of a number of recent studies, where 
the authors investigated its steady-state correlations [40-42] as well as its dynamical 
properties [43, 44]. 

In the rest of this section we derive the critical line of the canonical ABC model 
following the lines presented in [18]. We use, however, a different notation which is 
found to be convenient for analyzing the grand-canonical generalizations of the ABC 
model presented in the following sections. 

The weakly-asymmetric ABC model can be studied in the hydrodynamic limit 
where the probability to observe a particle of type a at site i is represented by a 
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continuous density profile, Pa{x), with x = i/L G [0,1) [38, 45, 46]. The canonical 
large deviation function of the hydrodynamic profile is given by 



Jo 



x) In Pa{x) + (1 — r) In (1 — r) 



+ /9 



dx 



/ dzy2iPa{x)pa-l{x + Z) - -1 

Jo ^ o 



(28) 



were p = {pA, Pb, Pc) and the index a runs cyclically over A,B and C. Equation 
(28) implies that the probability to observe a density profile, p{x), is given by 
P[p{x)] ~ e~^-^^P^^'>''^\ The first integral in (28) is the entropy of the density profile (the 
logarithm of the number of microstates corresponding to p), and the second integral is 
the continuum limit of the Hamiltonian (25). In order to find the critical point of the 
model we expand the density profile around the equal-densities homogeneous profile in 
terms of the following Fourier modes, 

/ 1 \ 



2ninx 



(29) 



where Oq = and ci^j = (a„,A, o,n,B, (in,c) is a complex vector which obeys = a* ^ 
in order to assure that p{x) is real. Here and below the ★ superscript denotes the 
complex conjugate. Since we study the model close to the homogenous phase we assume 
\0'n,a\ ^ 1- The lowest order terms in an,a the Landau expansion of (28) are given by 



Jf({a}) = r^r) 



(30) 



r 



i/3 
2Tm 



( 



\ 



1 
-1 

1 -1 



-l\ 



1 




i_„ + 0(a^_J, 



where the denotes the transposed vector of ci„, / denotes the identity matrix, and 
^ySj) = rln(r/3) + (1 — r) ln(l — r) is the free energy density of the homogenous phase 
multiplied by /5. 

The stability analysis of the Fourier modes can be carried out more conveniently by 
first diagonalizing the lowest order terms of J-'({a}). To this end we express the density 
profile in terms of the eigenvectors of the matrix which appears inside the square brackets 
in (30) given by 



u 





1 




1 


\ 






g-27ri/3 










g27ri/3 


/ 



1 

7! 



\ 



327ri/3 
-27ri/3 



The profile can thus be written as 

oo 

p(x) = ^u+ ^ [6„v + &*„v*]e 



(31) 



(32) 



where 6„ are the complex coefficients of the Fourier expansion. Since the steady-state 
profile of the inert vacancies is fiat, namely p^ix) = 1 — pa{x) — Pb{x) — pcix) = 1 — r is 
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Figure 4. Geometric representation of the amplitude bo which measures the deviation 
from equal densities. The angle in the sketch is given hy (p ~ arg(6o) +7r/6. The lengths 
of the dashed lines perpendicular to the edges denote the average densities of the three 
species, given by (33). 

independent of x, we omitted from (32) the amplitudes of the vector ue^'^*"^ for n 7^ 0. 
By integrating (32) over x, bo can be shown to set the deviation from equal-densities 
point through 

f f f Stt Stt 

{rA,rB,rc) = (^ , ^ , g) + |&o|(cos cos((^ + y ), cos(^ - y )), (33) 

where ip = arg(6o) +7r/6, as depicted in figure 4. By summing the vector components in 
(33) we see that 60 controls the differences between the densities, — and — rc, 
but not the overall particle density, which is given by r = ta + tb + tq- Since the 
canonical model undergoes a phase transition only for equal densities, we assume in this 
section that bo — 0. 

To lowest order in 6„ the canonical Landau expansion is given by 

+ 0{\bX). (34) 

At /3c = 2n\/3/r, the first mode, 61, becomes unstable, whereas the other modes are 
linearly stable. Minimizing the free energy with respect to bn close to the transition 
point, /3 — 27r-\/3/r <S 1, yields to lowest order 

6i+3Z bi-3i {b*f-'\ (35) 
for / = 1,2,... This scaling is a result of the third order coupling terms in J-' of the form 
^ni^n2^-ni-n2 ^ni^n2^-ni-n2 • '^^^ Other modcs Vanish, 

b±3i = 6_i±3i = 0, (36) 

since they are driven by 6_i which is not excited. The ABC model thus has a complex 
primary order parameter, bi, which drives some of the other Fourier modes. The analysis 
of Landau expansion of the ABC model thus requires a strait-forward generalization 
of the generic canonical analysis presented in section 2.1, where the primary order 
parameter, m, is real. By expressing the higher order Fourier modes in terms of bi and 



H{b}) = Mr) + E 



,3 



2'Kn 



)\bn\ 
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inserting the result back into (34) we obtain the Landau expansion of the model in terms 
1 61 1 alone, given by 

= Mr) + /2 W, r) + U iP, r) + O , (37) 

where 

The full derivation of /2 and /4 is found in [18]. We can now compare (37) with the 
generic Landau expansion, found in (3), and derive the phase diagram of the canonical 
ABC model following the hues presented for the generic case in section 2.1. Since /4 > 0, 
we conclude that the canonical model exhibits only a second order transition line defined 
by the equation f2 — 0, which yields 

= 27rV3/r. (39) 

In the following two sections we study the phase diagrams of two different grand- 
canonical generalizations of the ABC model. In section 3.2, Model 1 is shown to exhibit 
the same second order transition line as in (39). However, in contrast to the canonical 
ensemble, in Model 1 this line turns into first order at a tricritical point. In section 3.2 we 
study Model 2 which is shown to exhibit a second order transition line which is different 
than (39). Below we define the two grand-canonical generalization by performing the 
appropriate Legendre transform of T defined in (30). It is worth noting, however, that 
the two models can be studied by considering specific particle-nonconserving dynamical 
rules in addition to (23) and (24) that obey detailed balance with respect to the ABC 
Hamiltonian. These dynamical rules are discussed for Model 1 in [17, 18] and for Model 
2 m [19]. 

3.2. Model 1: Grand- canonical ABC model with ta — tb and tb — tc fixed 

In this section we consider a grand-canonical ABC model where the differences between 
the overall densities of the species, ta — tb and tb — rc, are fixed while the overall 
particle density, r — rA + tb + rc, is allowed to fiuctuate [17, 18]. We study the model 
for equal densities, = tb — rc, which is the only case where it exhibits a phase 
transition between the homogenous and ordered phases. The large deviation function 
of the density profile is given in this case by 

g [p{x) , r, //] = J- [p{x) , r] - /3pLr, (40) 

where /i is the chemical potential. We study the critical point of the model by expanding 
(40) close to the homogenous phase, while allowing for fiuctuations in r denoted by Sr. 
The density profile is thus given by 

p{x) = '-^u + J2 [&nv + 6^v^] e^— , (41) 
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where bo = since wc consider the equal-densities case. In terms of the generic model 
defined in section 2.1, it will be shown below that bi and 6r play the role of m and 5q , 
respectively. Using (41) we rewrite the Landau expansion obtained in [18] as 

em) « e.M + + (^ ^) 1^-1' - ^l".!'*^ 

+ (2(1^ + f*'-)^ - ^(*'*- + + JsIM' + O .(42) 

where Qh (^) = -^/i ~ /^A*^- The coupling term between the fluctuating field, Sr, and 
the primary order parameter, 61, in (42) is of the form 5r|6ip. Since this is compatible 
with the generic Landau expansion presented in (8), where f{m,q) — f{—m,q), we 
expect the canonical and grand-canonical ensembles to exhibit the same critical lines. 

The form of the coupling between bi and 6r can also be deduced from the symmetry 
properties of the model, without deriving (42) explicitly. The large deviation function, 
J^, is invariant under cyclic exchange of the species label, A — )■ S, 5 — > C and C — > A, 
given in terms of the Fourier modes by 

bn ^ ftne-^'^'/^ (43) 

In the equal-densities case J-" can be shown to be invariant also under a spatial shift 
of the density profile by a continuous parameter, A, which implies that J^[p{x),r] = 
T [p(modi(a; -f- A)), r]. In terms of the Fourier modes the translation of the profile takes 
the form of 

br, ^ fe^e-^'^^"^. (44) 

Since the overall particle density, r + 5r = Xla /o^ PqI^)*^^' invariant under (43) and 
(44), we expect the lowest order coupling between b\ and 5r in (42) to be of the form 
|6ip5r. 

Our conclusion that the grand-canonical critical line is identical to the canonical 
one can be demonstrated by calculating the former exphcitly. To this end, we first 
minimize Q with respect to and 6_2 in order to obtain their expression in terms of 
the primary order parameter, 61, 

{h\)\ br =^-^^^\b,\\ (45) 



r (4V37r + ^r) 

Inserting (45) into (42) yields the following Landau expansion 

Q{b^) = Gh (r) + 92 (/3, r) \b,\' + g, (/3, r) \b^\' + O {\b^\') , (46) 



where 



9 



r 27T 
V3(3r + 6n 3 (1 - r) 



V3/3r + 127r 



7^/4(Ar). (47) 
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Figure 5. Phase diagram of the canonical ABC model (a) and Model 1 (b), plotted 
in the (/i, T — 1//3) plane. The two diagrams show the same second order transition 
line (thin solid line). In the grand-canonical ensemble the transition line turns into a 
first order (thick solid line) at a tricritical point {-k). 



We can now compare (46) to the generic Landau expansion, found in (11), and derive 
the phase diagram of the grand-canonical ABC model following the hues presented in 
section 2.1. On the critical line, /3 — Pc — 27r-\/3/r, where g2 — 0, we find that 

94Pc,r) = ^-^^^^, (48) 

which is positive for r > 1/3. We thus conclude that the canonical and grand-canonical 
ensembles exhibit the same second order transition line at 

^c = 27rV3/r, for r > rrcp = 1/3. (49) 

Below the tricritical point (r < ttcp), where 5^4 < 0, we find a first order transition line 
in the grand-canonical ensemble, whereas the canonical transition line remains second 
order for all r. The resulting grand-canonical phase diagram is plotted in the (/.t, 1//3) 
plane in figure 5b. The figure shows a second order transition line (thin solid line) which 
turn into a first order line (thick solid line) at a tricritical point (^). The first order 
transition line in the figure is obtained from the exact expression of the steady-state 
density profile of the model, derived in [18]. For comparison we plot the canonical phase 
diagram in figure 5a, which is composed only of a second order transition line. The 
chemical potential in the canonical ensemble is given by /i — dJ-'h/dr. As expected, due 
to the symmetry properties of T, the ensemble inequivalence presented in figure 5 is 
consistent with the behavior of the generic model studied in section 2.1. 

3.3. Model 2: Grand- canonical ABC model with ta + tb + tq fixed 

In this section we analyze a different grand-canonical ABC model where r^, and rc 
are allowed to fluctuate independently, while maintaining the overall particle density, 
r = + tb rc, fixed. 



17 



The large deviation function of the density profile of Model 2 is given by 

G[p{x),ra,Xa\^ J^[p{x),r]- ^ XaVa: (50) 

a=A,B,C 

where Aq, are the conjugate fields which control the fluctuations of r^. Here we study 
case = As = Ac = /X, for which the model undergoes a phase transition. Since r is 
fixed, changing the value of A^ shifts ^ by a constant without affecting the resulting 
phase diagram. The field /i is thus an arbitrary parameter in this model, which can be 
set to zero. We expand Q in terms of small deviations from the homogenous profile, 
defined as 

oo 

pix) = + ^ [bnv + 6* „v*] e^-'^^ (51) 

n=— oo 

Unlike in the canonical model, where 60 is fixed, here it is allowed to fluctuate, leading 
to a significantly different phase diagram. To lowest order in 6„ we obtain 



Sm) = &(r) + (2 - |)|6„p + 



E 



,3 /3\/3,,, 2 /SVS 



' r 27rn 27m 



+ 0{bt). (52) 



Due to the bilinear coupling between 69 and 6„ in Q for arbitrary n, all modes have the 
same symmetry and any one of them can be viewed as the primary order parameter. As 
shown below, there is a specific combination of all the modes which becomes linearly 
unstable at the critical point. The bilinear couphng term suggests that this model is 
compatible with the generic Landau expansion in (15), where f{m,q) — f{—m,—q), 
with bi and bo playing the role of m and q, respectively. In Model 2, however, we find 
an additional infinite set of parameters, 6„ for n 7^ 0, 1, which are coupled linearly to 
60 and thus play a role in setting the critical point of the model [47]. Even though (52) 
has multiple order parameters coupled linearly to 69 and each of them is complex, the 
resulting behaviour should be similar to that of the generic model analyzed in section 2.2. 
This imphes that the canonical and grand-canonical ensembles should exhibit different 
critical points, as demonstrated explicitly below. 

The form of the coupling term between bg and bi in (52) can be deduced 
without performing the Landau expansion explicitly, simply by analyzing the symmetry 
properties of the model. Here, since we cannot restrict the system to remain in the 
equal-densities point, the translational symmetry is broken and J- is invariant only 
under cyclic exchanges of the species' labels given by 

bn ^ 6„e-2-/3_ (53) 

This implies that the Landau expansion of T may contain bilinear coupling terms of the 
form foo^'n ^^<i ^o^n- In principle, from these symmetry considerations, we also expect to 
find terms of the form bnb* with n,l ^ 0. However, explicit calculations show that the 
coefficients of these terms vanish. 
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We now demonstrate that the two ensembles exhibit different critical lines by 
computing the grand-canonical one explicitly. Since each of the Fourier modes can 
be regarded as the primary order parameter we choose to express them in terms of bo 
and then obtain an expansion of Q in powers of bo. This can be done by first minimizing 
Q with respect to 6„ for n ^ which yields 

-)bn - -7^ — ^^o = 0, (54) 



dbt r 2Tin 27m 



and thus 



bo. (55) 



/3r 

We insert (55) this into (52) and obtain a Landau expansion of Q in terms of bo, 



G{bo) = Gh{r) + 



,S /3v/3^ 1 



2 ^ 27r ^ Wl 2 _ ^ 



\bo\' + 0{\bo\'). (56) 



By simphfying the coefficient of |6oP above, as detailed in Appendix A, it can be shown 
to vanish at 

= 27rA/3(/ + \) I e Z. (57) 

The highest temperature, T = for which the instability sets in is thus obtained for 
/ = 0, yielding 

A=^. (58) 

At this point bo becomes unstable along with the rest of the modes which depend linearly 
on bo through (55). This critical point has been derived in [19] by computing the exact 
minimizing density profile of the large deviation function of Model 2 (50). 

As expected, since the Landau expansion involves a bilinear coupling term between 
the fluctuating field, bo, and the order parameter of the canonical model, bi, the canonical 
and grand-canonical models exhibit different critical points at /3 = 2-K\^/r and 
f3 = 27rv^/3r, respectively. The phase diagrams, derived within the two ensembles, are 
shown in figure 6. The figure corresponds to the two-equal densities case, rA = tb ^ tc, 
obtained by setting = = —Ac/2. In the canonical ensemble we find a second 
order phase transition point (•) on the A^ = line. Below this point the model 
undergoes a first order phase transition between a C-rich [ta — tb < tc) and a C- 
poor [ta — tb > Tc) ordered phases, depicted by the discontinuity in A^ = dT jdvA 
(vertical lines in figure 6a). In the grand- canonical phase diagram, shown in figure 6b, 
the second order phase transition (•) is found at a higher temperature. Below the critical 
temperature the model undergoes a first order phase transition at A^ = between a 
C-rich and a C-poor ordered phases (thick solid line). The C-rich (C-poor) phase is 
meta-stable for A^i > (A^ < 0) up to the dashed line where it becomes unstable. The 
stability limits (dashed lines) are obtained using the exact steady-state density profile 
of the ABC model, derived in [19, 38]. 
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Figure 6. Phase diagram of tlie canonical ABC model (a) and Model 2 (b), plotted 

in the {\a,T = plane for the case where = = —Ac/2 and r ~ 1. Here, 
the two ensembles exhibit different critical points (•). Below the critical points we 
find a first order transition, denoted in (a) by a discontinuity in A^ (vertical lines) and 
in (b) by a thick solid line. The dashed lines in (b) mark the stability Umits of the 
meta-stable phases, computed using the exact solution of the density profile. They 
define a coexistence region where a two different ordered profiles are locally stable in 
Model 2. See details in text. 



4. Conclusions 

In this paper we relate the order of the phase transition where two ensembles of a long- 
range interacting system may become inequivalent to the symmetry of the underlying 
order parameter of the transition. This is done within the framework of Landau theory, 
which constitutes a general approach for studying systems with long-range interactions. 

For a transition governed by an order parameter m and a field q which is allowed 
to fluctuate only in the higher ensemble, we consider two types of symmetries of the 
Landau free energy density, f{m,q). In the first case, where f{m,q) — f{—'m,q), the 
lowest order coupling between m and q is of the form m^q, yielding fluctuations in q that 
scale as (5g = g — g* ~ m^. These fluctuations cannot affect the critical line of the model 
which depend on the coefficient of in the Landau expansion of the model. We thus 
find that the two ensembles, where q is either fixed or allowed to fiuctuate, differ only 
in their first order transition lines and tricritical points which depend on higher order 
terms in the Landau expansion. In the second case, we consider a symmetry of the form 
f{m, q) = f{—m, —q) and analyze the stability of the phase where m* — q* — 0. Here, 
the Landau expansion may involve a bilinear coupling term of the form mq, yielding 
fluctuations in q that scale as 6q ~ m. These fluctuations affect the coefficient of the 
term in the Landau expansion, causing the two ensembles to exhibit different critical 
lines. 

An example for the first symmetry. /(m,g) = /(— m,g), can be found when 
considering m to be the average magnetization of the system and q to be the energy or 
the overall density, which are both invariant under m — )■ —m. Since this is the most 
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natural choice for q, our derivation explains why most of the examples of ensemble 
inequivalence that were studied in the past are of this type (see e.g. [13-18, 27]). The 
second symmetry, /(m, q) = f{—m, —q), is shown above to exist in an anisotropic mean- 
field XY model as well as in the generalized ABC model studied in [19]. We provide 
detailed examples of the two different symmetries, by deriving the Landau expansion of 
two distinct grand-canonical generalizations of the ABC model. Our general analysis is 
found to be applicable for the ABC model, even though the latter has a complex order 
parameter and a three-fold symmetry. 

We thus conclude that the symmetry of the model sets the order of transition 
where ensemble inequivalence may be observed. If the symmetry permits a bilinear 
coupling term in the Landau expansion between the order parameter which governs the 
transition in the lower ensemble, m, and the parameter which is allowed to fiuctuate 
in the higher ensemble, q, the two ensembles would generically exhibit different second 
order transition points. 
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Appendix A. Derivation of the critical point of Model 2 



We wish to find the critical point of Model 2, described in section 3.3. To this end we 
obtained the Landau expansion of the model in powers of bg, given by 

,3_^ /3^/3 
Y 2' 



E 



27r ^ Wl 2 



n 



\bo\' + 0{\bo\'). (A.l) 



The critical point is found where the coefficient of |6oP vanishes. This coefficient can be 
simplified by expressing it in terms of the DiGamma function, defined as 



(A.2) 



where 7 is the Euler-Mascheroni constant. Using the recursive properties of the 
DiGamma function given in [48] (section 6.3), 

i>{^-y) = i>{y)+^cot{7ry), (A.3) 

and 

(A.4) 



^'(l + Z/) = V'(?/) + -, 
y 

the sum over n in (A.l) can be written as 

1 



E 



- = -0(1 + ?/) --0(1 -y) = - -7rcot(7ry), 
n^/y — n y 



(A.5) 
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where y — ^r/27r-\/3. Equation (A.l) can thus be written as 



(A.6) 



V 2' 



y a/3 



-= I I TT cot Try 



|&or + 0(|6on. 



The critical point is found where tlie coefficient of |6oP vanishes, yielding the equation 



TT 



1 1 

TT cot Try = 

y y V3 



which can be written as 

I3r 



cot 



2VS 



1 

7^' 



(A.7) 



(A.8) 



(A.9) 



Equation (A.8) has several solutions given by 

f3r = 27iV3{l + -) / e Z. 
3 

The solution for / = has the highest temperature, T — 1//3, and thus corresponds to 
the critical point of the model, 

27rV3 



3r 



(A.IO) 



This point was derived in [19] by the studying the exact expression for the steady-state 
density profile of the model. Note that the critical point is different than the canonical 
one, given by /3c = ^^7^- They differ by a factor of 3, whose source is discussed in detail 
in [19]. 
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